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1 INTRODUCTION 



STATIC AXIALLY SYMMETRIC SOLUTIONS OF 
EINSTEIN- YANG-MILLS EQUATIONS WITH A NEGATIVE 
COSMOLOGICAL CONSTANT: THE REGULAR CASE 

Eugen Radu 

Albert-Ludwigs-Universitdt, Fakultat fur Physik, Hermann- Herder- Strafie 3 
Freiburg D-79104, Germany 

Abstract 

Numerical solutions of the Einstein- Yang-Mills equations with a negative cosmological constant are 
constructed. These axially symmetric solutions approach asymptotically the anti-de Sitter spacetime 
and are regular everywhere. They are characterized by the winding number n > 1, the mass and the 
non-Abelian magnetic charge. The main properties of the solutions and the differences with respect 
to the asymptotically fiat case are discussed. The existence of axially symmetric monopole and dyon 
solutions in fixed anti-de Sitter spacetime is also discussed. 

o 
O 
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After the discovery by Bartnik and McKinnon (BK) of a nontrivial particlelike solution of the Einstein- 
Yang-Mills (EYM) equations ED], there has been a great deal of numerical and analytical work on various 
aspects of EYM theory. A large number of self-gravitating structures with non-Abelian fields have been 
found (for a review see M). These include black holes with non-trivial hair, thereby leading to the 
possibility of evading the no-hair conjecture. 

Most of these investigations have been carried out on the assumption that spacetime is asymptotically 
flat. Less is known when the theory is modified to include a cosmological constant A which greatly changes 
the asymptotic structure of spacetime || . 

For a positive cosmological constant, the behavior of the solutions is similar in many respects to that 
of asymptotically flat geometries J4j . In particular the configurations are in both cases unstable || . 
q-( If we allow for a negative cosmological constant, the solution of the matter free Einstein equations 

possessing the maximal number of symmetries is the anti-de Sitter (AdS) spacetime. Being a maximally 
QQ, symmetric spacetime, it is an excellent model to investigate questions of principle related to the quantisa- 

tion of fields propagating on a curved background, the interaction with the gravitational field and issues 
related to the lack of global hyperbolicity. Also, lately has been a lot of interest in asymptotically AdS 
■ spacetimes, connected with string theory and related topics. 

Recently, some authors have discussed the properties of soliton and black hole solutions of the EYM 
system for A < (i.e. an asymptotically AdS spacetime |j, ^] ||). They obtained some surprising results, 
which are strikingly different from the BK type solutions. First, there is a continuum of solutions in terms 
of the adjustable shooting parameter that specifies the initial conditions at the origin or at the event 
horizon, rather then discrete points. The spectrum has a finite number of continuous branches. Secondly 
there are nontrivial solutions stable against spherically symmetric linear perturbations, corresponding 
to stable monopole and dyon configurations. The solutions are classified by non-Abelian electric and 
magnetic charges and the ADM mass. When the parameter A approaches zero, an already-existing branch 
of monopoles and dyon solutions collapses to a single point in the moduli space. At the same time new 
branches of solutions emerge. A fractal structure in the moduli space has been noticed 0]. 

As observed in JtJ, these solutions may have profound consequences in the evolution of the early 
universe. 

In Ref. 10, PH regular gravitating monopole and dyon solutions in Einstein- Yang-Mills-Higgs (EYMH) 



theory with Higgs field in the adjoint representation were shown to exist in asymptotically AdS spacetime. 
As happens in asymptotically flat space, a critical value for the Newton constant exists above which no 
regular solution can be found. The presence of a cosmological constant enhances this effect, the critical 
value being smaller than the value found for A = 0. 
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The global existence of a solution of Cauchy problem for the YMH equations in AdS spacetime is 
discussed in |l2| . 

However, so far only spherically symmetric solutions have been found. 

For A = 0, a SU(2) YM theory coupled to a scalar Higgs field |]| |]|, dilaton |]| or gravity |1| is 
known to possess also axially symmetric finite energy solutions. In the 1990s, the numerical calculation of 
these configurations was one of the most important developments in the domain. 

A natural question arises: do non-Abelian axially symmetric solutions also exist for a nonzero cosmo- 
logical constant? And if this is the case, how does the nonasymptotically flat structure of spacetime affects 
these configurations? 

The aim of this paper is to address the above questions for the SU(2) gauge group, under the assumption 
of axial symmetry, for an asymptotically AdS geometry. 

The corresponding problem for a vanishing cosmological constant and a purely magnetic gauge field 
has been exhaustively discussed in |L6|. Representing generalizations of the BK solutions the solutions 
obtained by Kleihaus and Kunz in |lq ] have no nonabelian charges but are characterized by two integers. 
These are the node number k of the gauge field functions and the winding number n with respect to the 
azimuthal angle (p. The spherically symmetric BK solutions have winding number n = 1. A winding 
number n > 1 leads to axially symmetric solutions. As discussed in Jl6[ | , these regular axially symmetric 
solutions have a torus-like shape. With the z— axis (0 — 0) as symmetry axis, the energy density has a 
strong peak along the p axis {9 = 7r/2) and decreases monotonically along the z— axis. 

The solutions we are looking for are the asymptotically AdS analogues of these configurations. 

Although some common features are present, the results we find are rather different from those valid 
in the A = case. A different behavior is noticeable especially for the lower branch of axially symmetric 
solutions. These distinctions arise from differences that already exist in the spherically symmetric case. 

A negative A implies a continuum of axially symmetric solutions. For a fixed winding number, they can 
be classified in a finite number of branches and have continuous values of mass and non-Abelian magnetic 
charge. The radial and angular dependence of the metric and gauge functions can also be different from 
the case discussed in 

The paper is structured as follows: in the next section we explain the model and derive the basic 
equations, while in Section 3 we discuss solutions of YM equations in a fixed AdS background. These 
solutions were found to be important when discussing the scaling properties of the mass spectrum ||] and 
have no flat space couterparts. Some features of the axially symmetric solutions possessing a net YM 
electric charge in a fixed AdS background are also presented in this Section. The general properties of 
the axially symmetric gravitating solutions are presented in Section 4 where we show results obtained by 
numerical calculations. We give our conclusions and remarks in the final section. 

2 GENERAL FRAMEWORK AND EQUATIONS OF MO- 
TION 

2.1 Einstein- Yang-Mills action 

The basic equations for a static, axially symmetric SU(2) gauge field coupled to Einstein gravity (without 
cosmological term) are well known (for details see e.g. (l6|). Here we derive them for a nonzero A, without 
going into details. We will follow most of the conventions and notations used by Kleihaus and Kunz (KK) 
in their papers. 

The starting point is the Einstein- Yang-Mills action 

S = J d^xV^ij^ (n - 2 A) - irr( V")] , (i) 

where the field strength tensor is 

Ffj, v = d^Au - d v A„ + ie [A^jA^,] , (2) 
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and the gauge field 

A, = \r a Al. (3) 
Variation of the action (fil) with respect to the metric g^ u leads to the Einstein equations 



R)j,v - t:9 + kg^v = SnGTuv, (4) 



where the YM stress-energy tensor is 



Tfj, u = 2Tr(F^^ 5 ^ - \ 9llv F a0 F^). (5) 

Variation with respect to the gauge field leads to the matter field equations 

V^" +ie[A li ,F^\ =0. (6) 

2.2 Static axially symmetric ansatz and gauge condition 

We generalize the isotropic axisymmetric line element considered by Kleihaus and Kunz in |l6)| for a 
nonzero A 

ds 2 = -/(l - ^)d? + ^(-^r + rW) + ~r 2 sin 2 9d<t>\ (7) 
3 / 1 - J 

where the metric functions /, to and 2 are only functions of the coordinates r and 9. 

A suitable parametrization of a purely magnetic Yang-Mills connection in terms of spherical coordinates 

is 

A r = 2^iM)r£, 

Ae= ±(l-H 2 (r,6))Tl 
2e v 



As = 



T '<- rr-'l' 

H 3 (r 7 9)^- + (l-H 4 (r,8))^- 



(8) 



Here the symbols r™, t@ and denote the dot products of the cartesian vector of Pauli matrices, r = 
(Ti,T2,Ta), with the spatial unit vectors 

e™ = (sin 9 cos n</>, sin sin n<j>, cos 6*) , 
e 'g = (cos cos n<j), cos sin n0, — sin 9) , 

e£ = (— sin n(j), cos n<p, 0), (9) 

respectively. This ansatz is axially symmetric in the sense that a rotation around the z— axis can be 
compensated by a gauge rotation. It satisfies also some additional discrete symmetries [ijj]. To fix the 
residual abelian gauge invariance we choose the usual gauge condition |l6f| 

rd r Hi - d e H 2 = 0. (10) 

This ansatz contains an integer n, representing the winding number with respect to the azimutal angle (p. 
While if covers the trigonometric circle once, the fields wind n times around. Note that the spherically 
symmetric ansatz corresponds to n = 1, Hi = H 3 = 0, H 2 = H 4 = w(r). 
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2.3 Field equations 

From (f§ and (§ we obtain a set of seven nonlinear elliptical partial differential equations which can be 
solved numerically. 

Within this specific ansatz and the gauge condition (QOh we derive the matter equations 



r 2 Hi, r , r + Hififi + H 2 ,e - n 2 j (rH 4 , r H 3 - rH 3 . r H 4 + (iff + H\ - l) H x ) J sin 2 
rH 2>r + H he -n 2 — (2H 4 H 3 + rH A , r )J sin6»cos6» - n 2 —H 4 



' fNfl \ 



sin 2 6[Hie + rH 2 r ] In — + sin 2 6 [rH 1 r - H 2 e] r In 



r 2 H 2 ^ r + H 2>e ,e - H he ) +n 2 ^ {H 3 H 4 , e - H 4 H 3fi - (iff + H 2 - 1 ) H 2 ) J sin 2 
H 2fi - rH l>r + n 2 — (-2H 2 H 3 + H 4j6 )J sin cos + n 2 — {H 4 - H 2 ) 



ffNVl\ 

\ m J.r 



(11) 



l^yl) +sin 2 9\H i,e+ rH 2 . r lrLa f 



f 

+ j^H 3 . e , 


6 — 


- j sin 2 9 H 


1 

' N 



sin 2 6 [rHx, r - H 2fi )] In — + sin 2 6 [Hi, 9 + rH 2 , r ] r In 



2 7,, , + -II ;,<.,< II, ( Hi + ^ + H 4 (H 4 - 2rff 4 ,r) - #4 (r£Ti, P - ^ 

2i/ ^ j sin 2 + 1 (ff 8 ,9 + M 4 - A^ 2 - sin cos - ^ 
■i sin 2 [iJ 3 ,e - 1 + H 2 H 4 + cot 9H 3 ] In ^ j + sin 2 6 [rH 3 , r - HiH A ] r In (J^=- 



(12) 



(13) 



r'H^r + -H 4 . efi - H 4 l^Hf + J ~ H ^ ~ (-1 + 2H 3,o) + 2rH 3 , r H 4 

-H 3 (rH ltr - ^j^j sin 2 6 + (-Jf 2 /i 3 - #1^ + rfTi, r JV - ff 2 ,e + #4,6>) sin (9 cos 6» 
H 2 -H 4 , f _ ; . 2/irrr rr rr ,„ „,__ i/ni _/'/ 



sin 2 [riJ 4 , r + #i#3 + cot 0iTi] r In ( ^ ) , (14) 



y , -sin^[JJ 4 ,e-if3^2-(^2-^4)cot0]lnl ^ 
where iV = 1 — -|-r 2 . The resulting equations for the metric functions are 



\{ NrJ f r -T + I -f l i))- ! Y^f + 2 ^) + 2A( 7 " 1)1 (15 » 
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_ _ 1 ( /' Nr 2 m^ r + Nrm ir + m t e t e /ra,e\ 2 / rm^ \ 2 \ Jfi; 

J V r ' * + J ~ 4r 2 I ^ to V m ) V to ) J I 



-2 |'^V _ (Li) 2 + 2N (^ L + r -^A + N-™ r H ' r m ' 6 l ' 8 



f J \ I J \ I m J ml ml 

-^+2A(=-l). (17) 
In the above relations, the components of the energy-momentum tensor are 

j<r 1 ffY N F * | 1 F 2 1 / 2 p2 



2 VW r 2 rW 2 mZr 2 sin 2 r * 2 mZr 4 sin 2 



Tfl 2 (, m J r 2 Fre 2 rmY 2 sin 2 Fr + + 2 mZr 4 sin 2 9^ 

T* = KIXHf* i 1 PN F 2 i 1 /2 F 2 
* 2W r 2 r ^ 2 to/t 2 sin 2 r0 2 rmY 4 sin 2 9 H 

rpt = 1 ffY N F 2 1 / 2jV F 2 ■ 1 P F 2 

1 2\m r 2 r9+ 2 mlr 2 sin 2 r * 2 mZr 4 sin 2 



where, similar to |16| we define 

*% = 4 + r #2,r) 2 + ( ri/ i^ - ^) 2 ] . 



o . 2 

n sin I 



FA = n^sin^ 



[(rff 3 , r - i?ii? 4 ) 2 + {rH A , r + H t H 3 + cot 9H 1 ) 2 ] 



{Hi.e - H 2 H 3 - cot 9(H 2 - H±) f + (H 3 ,e - 1 + H 2 H 4 + cot 9H 3 ) 



2.4 Boundary conditions 

To obtain asymptotically AdS regular solutions with finite energy density the metric functions have to 
satisfy the boundary conditions 

d r f\ r =o = d r m\ r=0 = d r l\ r=0 = 0, (18) 

f\r=oo = Tll\r— oo = l\r— oo — 1- (1^) 

We assume also the flatness condition |18| 

m\e=o = l\e=o- (20) 
The boundary conditions satisfied by the matter functions are 

H 2 \r=Q = H^\ r=Q = 1, i?l| r= o = #3^=0 = 0, (21) 

at the origin and 

H 2 \ r=OQ = i?4|r=oo = Wo, fll|r=oo = #3|r=oo = 0, (22) 
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at infinity, where wq is an undetermined constant. For a solution with parity reflection symmetry, the 
boundary conditions along the axes are 

Hl\g = Q i7r /2 = H 3 \ e=Qn / 2 — 0, 

deH 2 \0 = o t7r /2 = d g H i \ e=0 T T /2 — 0, 

def\e=o,7r/2 = dem\e= ,7r/2 = d e l\ g=0: n/2 = 0. (24) 
Therefore we need to consider the solutions only in the region < 9 < 5. Regularity on the z— axis 
requires 

-H2|fl=o = #4|e=o- (25) 
Dimcnsionless partial differential equations arc obtained by the following rcscaling r — ► (V4ttG / e)r, A — * 
(e 2 /47rG)A. In Bjoraker-Hosotani conventions, this corresponds to taking a unit value for the parameter 
v = AnG/e 2 §. 

2.5 The mass and the Yang-Mills charges of the solution 

At spatial infinity, the line element (0) can be written as 

ds 2 = ds% + h^dx**" , (26) 

where are deviations from the background AdS metric cLsg. Similar to the asymptotically flat case, one 
expects the values of conserved quantities to be encoded in the functions h^ v . The construction of these 
quantities for an asymptotically AdS spacetime was addressed for the first time in the eighties (see for 
instance Ref. 19, 2Q|). However, the generalization of Komar's formula in this case is not straightforward 



j a = — ( P ^[G^tetViffw - huViti) + 2^4], (27) 



and requires the further subtraction of a background configuration in order to render a finite result. 

Using the Hamiltonian formalism, Henneaux and Teitelboim p0[ have computed the mass of an asymp- 
totically AdS spacetime in the following way. They showed that the Hamiltonian must be supplemented by 
surface terms in order to be consistent with the equations of motion. These surface terms yield conserved 
charges associated with the Killing vectors of an asymptotic AdS geometry. The general expression of a 
conserved quantity is 

1 

16tt 

where G ijkl = i(-5) 1/2 (3 <f +9 < 9 - 2.9 'g ), ^ is the component of the Killing vector ( A in the 
direction of the unit normal to the hypersurface t = const., hik is the deviation from the AdS metric, ~K % k is 

o 

the canonical momentum and is the covariant derivative with respect tot the background three-metric 

o 

9ik- The total energy is the charge associated with the Killing vector d/dt. 

It has been also checked that the Henneaux- Teitelboim energy, the Brown- York energy |2lJ and the 
Einstein and Landau-Lifshitz energy (l9) all agree for the Kerr-AdS spacetime. However, for the same 
spacetime, the generalized Komar mass has a different value f22|| . 

In this work we use the Henneaux- Teitelboim formalism to compute numerically the mass-energy of a 
gravitating EYM configuration. Substituting the ansatz (fjj) in ( |27| ) yields the following expression for the 
total mass 

M = lim A r /2 de ^0 (-^1 + '(Jtt+i)) r\ (28) 
r^oo 12 J \ rf dr f J 



The dimensionless mass is obtained by using the rescaling M — > (eG / '\/AkG)M . 

Solutions of the field equations are also classified by the nonabelian electric and magnetic charges Qe 
and Qm- The nonabelian charges defined by 

Qe - ' f dS k V=9Tr fffeoVr (29) 



are conserved because the Gauss flux theorem 
Within our ansatz Qm = n(l — to 2 ), Qe = 0. 
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3 YANG-MILLS FIELDS IN FIXED ADS BACKGROUND 



Because the asymptotic structure of geometry is different, in an AdS spacetime one does not have to couple 
the YM system to scalar fields or gravity in order to obtain finite energy solutions. Here the cosmological 
constant breaks the scale invariance of pure YM theory to give finite energy solutions. 

It is the purpose of this section to present both analytical and numerical arguments for the existence 
of nontrivial monopoles and dyons solutions of pure YM equations in a four-dimensional AdS spacetime, 
gravity being regarded as a fixed field (i.e. f = I = m = 1; also, in this section we don't use the above 
discussed rescaling). 

Although being extremely simple, nevertheless this model appears to contain all the essential features 
of the Bjoraker-Hosotani solutions. In this way, what might seem surprising (e.g. the existence of stable 
solutions) finds a natural explanation. 

The existence of these nongravitating solutions has recently been noticed in Q when discussing the 
scaling behavior of the EYM monopoles and dyons. 



3.1 Spherically symmetric solutions 

We start by briefly discussing the solutions obtained for n = 1 within the ansatz (||). In this case, the 
equation of motion has the simple form 

(^'(l-^ 2 ))' = 1 r2 i (30) 

where the prime denotes derivative with respect to r. The numerical results show the existence of a 
one-parameter family of solutions regular at r = with the behavior familiar from the gravitating case 

w(r) = l-br 2 + 0(r 4 ), (31) 

where b is an arbitrary constant. The asymptotic expansion at large r is 

w = w + w i - -I , (32) 

r 

where wq, wi are constants to be determined by numerical calculations. These boundary conditions permit 
a non- vanishing magnetic charge Qm- 

The overall picture we find is rather similar to the one described in H where gravity is taken into 
account. By varying the parameter b, a continuum of monopole solutions is obtained. As a new feature, 
we notice the existence of zero- and one-node monopole solutions (k = 0, 1) only. Also, for a fixed value 
of A, we obtain finite energy solutions for only one interval in parameter space, b m i n < b < b max (with 
bmin < 0). The energy of the solutions is an increasing function of the absolut value of b and diverges at the 
extremities of the interval. The allowed values of b corresponds approximately to the lower branch found 
in [Q when coupling to gravity. This is not an unexpected feature. We recall that, given a flat spacetime 
soliton solution, one can expect that it will have (asymptotically flat) gravitating generalizations. Apart 



from the fundamental gravitating solutions a sequence of radial excitations is likely to exist [j2|, 23 . For 
example, in flat spacetime the YMH system (with a doublet scalar field) has only one-node solutions; when 
gravity is included, the solutions exist for all k . 

For large enough r, it is possible to write an approximate expression of w(r) in terms of elliptic functions. 
A nontrivial exact solution of the YM equations is 

u= l/(l-A/3r 2 ) 1/2 , (33) 



describing a monopole in AdS spacetime with unit magnetic charge and mass VH5A)7r/8e 2 . It should be 
mentioned that this is not a new result. The existence of this exact solution has been noticed for the first 
time in pq| for a positive cosmological constant and a different coordinate system. 
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The arguments presented in Q for the (linear) stability of the nodeless n = 1 monopole solutions apply 
directly to the nongravitating case. In Bjoraker-Hosotani analysis (Ref. 0, Section VII), this corresponds 
to taking v = 0, H = 1 — Ar 2 /3, p = 1, while i5iJ = <5p = <5m = 0. We starts with the most general 
expression of a spherically symmetric YM connection 

J ^ = 2^{ n ( r ' t) T 3dt + v(r, t)r 3 dr + {w(r, t)ri + w(r, t)T 2 )dO 

+ (cot6»T 3 + w(r, t)r 2 - u)(r, t)n) sin^j. (34) 

All field variables are written as the sum of the static equilibrium solution whose stability we are investi- 
gating and a time dependent perturbation. In examining time-dependent fluctuations around finite energy 
solutions it is convenient to work in the Su(r, t) — gauge. By following the standard methods we derive 
linearized equations for Sw(r,t), Sw(r,t) and <5^(r, t) . The linearized equations for Sw(r,t) and <5z/(r, t) are 
obtained in for the more general gravitating case (as usual 5w(r, t) is determined by 5v[r, t)). For a 
harmonic time dependence e lM , the linearized system implies two standard Schrodinger equations. The 
analysis of the potential's properties in these Schrodinger equations can be done following Ref. 0. The 
standard arguments presented by Bjoraker and Hosotani are still valid and imply that for nodeless solutions 
there are no negative eigenvalues for Q 2 and thus no unstable modes. 



3.2 Axially symmetric solutions 

In addition to these spherically symmetric solutions we study their axially symmetric generalizations. 
Subject to the boundary conditions jl^)-(p^), we solve the YM equations numerically. 

Our methods are similar to those used by the authors of jl6| in their works. The KK scheme solves 
the field equations following an iteration procedure. One starts with a known spherically symmetric 
configuration and increases the winding number in small steps. The field equations are first discretized on 
a noncquidistant grid and the resulting system is solved iteratively until convergence is achieved. In this 
scheme, a new radial variable is introduced which maps the semi-infinite region [0, oo) to the closed region 
[0, 1]. Thus the region of integration is not truncated and the model converges to a higher accuracy. There 
are various possibilities for this transformation, but a choice which is flexible enough was x — ^r^, where c 
is a properly chosen constant. Typical grids have sizes 170 x 30, covering the integration region < x < 1 
and < 6 < tt/2. 

The numerical calculations are performed by using the program FIDISOL, based on the iterative 
Newton- Raphson method. Detailes on the FIDISOL code are presented in |^| . To obtain axially symmetric 
solutions, we start with the n = 1 solution discussed above as initial guess and increase the value of n 
slowly. The iterations converge, and repeating the procedure one obtains in this way solutions for arbitrary 
n. The physical values of n are integers. The numerical error for the functions is estimated to be lower 
than 10~ 3 . 

The energy density of these nongravitating solutions is given by the it-component of the energy mo- 
mentum tensor T"; integration over all space yields their mass-energy 

e = J {{/•;;/;; + v^ 3 * = ^ J\ 2 dr J\medel^F^ + \f*f^. (35) 

The energy and the magnetic charge of solutions are proportional to their winding number n. 

For every spherically symmetric configuration (i.e. a given value of parameter 6) we have obtained 
higher winding number generalizations. Moreover, the branch structure noticed for n = 1 is retained for 
higher winding number solutions. These solutions have very similar properties with the corresponding 
EYM counterparts. Therefore, the general picture we present here applies also in the next section. 

In Fig. 1 we present the gauge functions Hi and the energy density e as function of the radial coordinate 
r for the angles 9 = 0, tt/A and 7r/2 for three different solutions. Here the winding number is n = 3 and 
A = —0.01; also the mass is given in units Air/e 2 . The configuration with (Qm = 3, M = 0.332) represents 
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a higher winding generalization of the exact solution (p3|)(z.e. b — 0.00166 and wq — 0); we suspect the 
existence of a general analytic form of this solution (valid for n > 1). 

The configurations with ojq ^ have been arbitrarily selected; the solution with (Qm — 2.391, M = 
0.576) is obtained starting with a spherically symmetric solution with b = 0.003, while for (Qm = —13.284, 
M = 1.118) we have b = -0.001. 

The functions H 2 and H4 have a small 9 dependence, although the angular dependence of matter 
functions generally increases with Qm- We notice also that the gauge field function Hi remains nodeless 
and for every solution with wq > 1 it takes only negative values (Hi and H 3 are zero on the axes in Fig. la 
and lc). 

In Fig. 2 we show the energy density e and the gauge functions Hi for a nodeless solution with n = 3, 
k = 0, total mass M = 5.481 (in units Air/e 2 ), magnetic charge Qm = —37.185 as a function of the 
compactified coordinates p — x sin 9 and z — xcos9 (for A = —0.01) (here x = r/(50 + r)). The parameter 
b for the corresponding spherically symmetric solution is b = —0.0015. In this case, the functions Hi and 
H4 are almost spherical. The function Hi does not possess a non-trivial node and takes only negative 
values. 



3.3 Nongravitating dyon solutions 

The existence of dyon solutions without a Higgs field is a new feature for AdS spacetime M . If A > the 
electric part of the gauge fields is forbidden 0, |27j. In order for the boundary conditions at infinity to 
permit the electric fields and maintain a finite ADM mass we have to add scalar fields to the theory. 
The YM ansatz (JsJ) can be generalized to include an electric part (see e.g. G3) 

A t = H 5 (r,e)^+H 6 (r,9)^-. (36) 
2e 2e 

A possible set of boundary conditions for the electric potentials H§ , H§ is 

H 5 \ r=0 = H e \ r=0 — 0, H 5 \ r=OQ — u , H$\ r=oc — 0, (37) 

and 

deH 5 \ g=Q ^/2 — H 6 \ g=0t7r /2 = (38) 
for a solution with parity reflection symmetry. The equations of motion in this case are 

= (^r 2 H hr , r + Hi fifi + H 2fi - n 2 (rH^ r H 3 - rH 3 , r H 4 + (iff + H 2 - l) Hi) ^ sin 2 9 
+ (rH 2 , r + H h9 - n 2 (2H X H 3 + rH i>r )) sin0cos0- n 2 Hi 

+ sm 2 9[rH hr -H 2 , e }r\n(N) r + - [rH 5>r H 6 - rH 6 , r H 5 + H X (H 2 + g 6 2 )] , 

= (^r 2 H 2 ,„. + H 2fi>e - H lfi ) + ^ (H 3 H it8 - H^.g - (H 2 + H 2 - l) H 2 ) ^ sin 2 9 
H 2 . e - rHi. r + — (-2H 2 H 3 + H ifi ) J sin 9 cos 9 + — (H A - H 2 ) 



sin 2 9r 2 



+ sin 2 9 [H he + rH 2>r ] r In (N) >r + [H 5 H 6 , e - H 6 H 5 , e + H 2 (H 2 5 + tf 2 )] , 



r 2 H 3 , r , r + ^H 3<e ,e - H 3 (hI + ^\ + Hi (H 4 - 2rH^ r ) - H A (rH hr - ^ 



(39) 



(40) 
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sm2 6 + jf + H ^ NH " H ") sin 9cos9 if 

sin 2 9r 2 H 

+ sin 2 [rfr a>r - HxH 4 ] r In (tf) iP + ^— £ [ff 4 ff 5 + #6 (#3 + cot 0)] , 



(41) 



-fls ( rJT ljr - ^) ) sin 2 9 + 1 (-# 2 F 3 - ^#1 + rffi, r JV - iT 2 , e + ff 4 , e ) sinflcosf ''" " 



N I / N ° 1 1 1 1 N 



sin 2 9 [rHi.r + HxH 3 + cot ftffj r In (7V) r + -5— £ [ff 4 ff 5 + ff 6 (ff 3 + cot 6»)] , 



(42) 



r 2 



+H 5 (NH 2 + Hi + n 2 H 2 )"jj sin 2 9 + A. (#5,0 - #5(^2 + « 2 # 4 )) sintfcostf, 

(43) 

= ^r 2 iT 6 , r , r + M - rF 1)r iT 5 - 2rH^ r H x + A (2H 6 , 6 H 2 + H 5 {H 2fi - n 2 H 3 H±) - H 6 (NH 2 

+H 2 + n 2 H 2 - 1)) sin 2 9 - + j- (h 6 , 6 + H 2 H 5 + n 2 (2H 3 H 6 + H 4 H 5 )^j sin 9 cos 9. (44) 

Similar to the asymptotically flat case [p8f , a vanishing uo implies a purely magnetic solution. To prove 
this, we express the electric part in (BF 



Ee = \ J F? t F? t ^-gd*x (45) 

as a surface integral at infinity. We use also the existence of the Killing vector d/dt which implies 

Fa = DiA t (46) 
and the YM equations (0). Thus we obtain the general result 

-E e = Tr(J{D t (A t F H ^)-A t D l (F lt V=g)d 3 x}) (47) 

and, for a regular configuration 

E e = Tr((f A t F rt dS r ). (48) 

Therefore, for our ansatz 

E e = ^ lim / r 2 d r H 5 sm9d9. (49) 
This result provides also an useful test to verify the accuracy of the numerical calculations. 
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When taking n — 1, H§ — u(r) and Hq = we find spherically symmetric nongravitating dyon 
solutions. In this case we have 



u(r) 




(50) 



at the origin (where a and b are arbitrary constants) and 



1 

u = uo + ui~ H 

r 



(51) 



at large r, where uq, tti, u>o, u>i are constants to be determined by numerical calculations. The expansion 
for w(r) is still valid. These boundary conditions permit non-vanishing charges Qm and Qe- In order 
to obtain the value of these charges at some distance, we have calculated the integrand in (^9|) in the 
numerical code. 

If the shooting parameter a is nonzero, we find dyon solutions. Solutions are found for a continuous 
set of parameters a and 6; for some limiting values of these parameters solutions blow up. Given (a, b) 
the general behavior of the gauge functions w, u is similar to the gravitating case; there are also solutions 
with Qm = but Qe =/= 0. The surprising numerically properties noticed for the gravitating case |Q are 
found also for our solutions (for example, Qm — — l/v^r at b = 0.0061 independent of the value of a). In 
fact we suspect that these properties reside in the nongravitating sector of the theory. 

When studying dyon solutions we notice the existence of higher node (k > 1) configurations. Also, 
there are solutions where w does not cross the r axis. For a fixed value of 6, the number of nodes is 
determined by the value of the parameter a. Typical spherically symmetric solutions are displayed in Fig. 



The same numerical method described above is used to obtain higher winding number dyon solutions. 

Axially symmetric dyon solutions are also known to exist in a SU(2) Yang-Mills- Higgs (when working in 
flat spacetime background). In the Prasad-Sommcrfcld limit, they are known analytically, while for finite 
Higgs self-coupling they have been recently constructed numerically jL4). To the author's knowledge there 
are no known regular axially symmetric dyon solutions (analytical or numerical) in a nonflat geometry. The 
toroidal shape of the energy density of the monopole solutions is retained for the n > 1 dyon solutions, as is 
illustrated in Fig. 4g. As a typical axially symmetrical configuration, we show in this three-dimensional plot 
the gauge functions Hi and the energy density e for the a solution with n = 2, k = 1, total energy E — 1.046 
(in units 47r/e 2 ), nonabelian charges Qm = —4.124 and Qe — —4.362 as a function of the compactified 
coordinates p = x sin# and z — x cos 6 (for a better visualization we define here x = r / (100 + 7')). The value 
of the cosmological constant is A = —0.01. This solution has been obtained starting from a spherically 
symmetric configurations with the shooting parameters b = 0.01 and a = 0.005. As seen in Fig. 4a-f, the 
gauge functions H2, -H4, H§ do not exhibit a strong angular dependence. The Hi, H3 functions remain 
nodeless while the second electric potential ifg always presents a complicated nodal structure and angular 
dependence. 

The effect of the presence of the YM electric charge is seen in Fig. 5, where the energy density is 
shown as a function of the radial coordinate for several values of the angle 9, both for an axially symmetric 
dyon solution and a monopole solution with the same values of magnetic charge and winding number 



4 AXIALLY SYMMETRIC SOLUTIONS IN THE PRESENCE 
OF GRAVITY 



In their paper |7j Bjoraker and Hosotani have used Schwarzschild-likc coordinates with a line element 



3. 



(A 



0.01). 




(52) 
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where 

= (53) 
r S 

Since we want to use the spherically symmetric solution as the starting point for the calculation of axially 
symmetric configurations, it is appropriate to transform the known Bjoraker-Hosotani solution to the 
coordinates which appear in our general metric (Q). 

Given the presence of a cosmological constant, this coordinate transformation is more complicated then 
the transformation in [fL6l. 



4.1 Coordinate transformation 

By requiring / = m and the metric functions / and m to be only functions of the coordinate r, the axially 
symmetric isotropic metric @j reduces to the form 

ds 2 = -/(l - ^r 2 )* 2 + ^(-^ + r 2 (dO* + sin 2 OA?)). (54) 



The relations ©-(M yields 



and 



f(r) _ r 



2 



- ~ = ft 1 (55) 



m(r) f 2 
dr 1 dr 



(56) 



Since the mass function m(f) is only known numerically, we have to numerically integrate ( pq ) to obtain 
r(f). Therefore we find 

m = TT | i¥ ,, (57) 

with 

r f°° i ( i 

^(f) = exp 



1 df' 



r \VH ) 

The integration constant is adjusted such that at infinity ft = 1, i.e. r = r. The integrand in ([58|) is well 
behaved at the origin, since m(r) ~ f 3 0. 

For a first branch solution, the values of the metric functions are close to one. Fig. 6a demonstrates the 
coordinate transformation for higher branch spherically symmetric solutions with k = 1 — 3 and different 
magnetic charges. The metric functions /, m and the gauge field function w are shown in figure Fig. 6b. 
These solutions resemble those obtained for A = 0, with quantitative differences only. For example a 
smaller value of metric function m at the origin has to be noticed. 



4.2 Numerical method 

We employ the same numerical algorithm as for the YM solutions in fixed AdS background presented 
above. To obtain axially symmetric solutions, we start with a n = 1 EYM solution as initial guess and 
increase the value of n slowly (for a fixed u)q). A second procedure, also employed for the first branch 
solutions, is to start with a known axially symmetric YM solution (with n = 2, 3, ...) as an initial guess for 
the full system. 

The numerical error for the functions is estimated to be on the order of 10 -3 or lower for first branch 
solutions and 1CP 2 in rest. This error depends also on the magnetic charge and mass of the solutions. 
Axially symmetric generalizations of the n = 1 solutions with a large mass are difficult to obtain. 

A set of A = test runs was carried out, primarily designed to evaluate the code's ability to reproduce 
the KK results. In this case, we have obtained an excellent agreement with the results of 
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4.3 Properties of the solutions 

For all the solutions we present we take a cosmological constant A = —0.01, also the main value considered 
in Ref. ffi]. However, a similar general behavior has been found for other negative values of A. 

Starting from a spherically symmetric configuration we obtain higher winding number generalizations 
with many similar properties. For a fixed winding number, the solutions can also be indexed in a finite 
number of branches classified by the mass and the non-Abelian magnetic charge. This is in sharp contrast 
to the A = case, where only a discrete set of solutions is found Jl6| . 

The metric functions f,m,l are completely regular and show no sign of an apparent horizon. 

We begin with a description of the lowest branch axially symmetric regular solutions. In this case the 
winding number is n > 1 and nodeless or one-node solutions are allowed. These solutions are of particular 
interest because they are likely to be stable against linear perturbations (for k = Q). As expected, the 
gauge functions Hi looks very similar to those of the corresponding (pure-) YM solutions. The general 
picture presented in Fig. 2 is valid in this case too. The typical values of the metric functions m, /, I are 
closed to one. These functions do not exhibit a strong angular dependence, while m and I have a rather 
similar shape. 

To see the change of the functions for an increasing n, we exhibit in Fig. 7 first branch solutions with 
k = 1, ujq = —3, and n = 1, 2 and 3. In Figs. 7a-d the gauge field functions are shown, in Figs. 7e-g 
the metric functions, and in Fig. 7h the energy density of the matter fields. These two-dimensional plots 
exhibit the r dependence for three fixed angles 8 = 0, ir/4 and ir/2. Note that the Hi, H$ functions remain 
nodeless (Hi and H3 are zero on the axes in Figs. 7a, c as required by the boundary conditions (p3|)). 
As expected, the angular dependence of the metric and matter functions increases with n. However, the 
location of the nodes of the gauge field functions H2, Hi does not move further outward with increasing 
n. At the same time the peak of the energy density along the r-axis slightly shifts outward with increasing 
n and increases in height. At the origin the values of the metric functions decreases with n. 

This behavior contrasts with the picture obtained in an asymptotically flat spacctime. 

In Fig. 8 the mass of the first branch solutions M is plotted as a function of the nonabelian magnetic 
charge Qm for various winding numbers. For the studied configurations we find that the total mass of 
a gravitating solution has a smaller value than the mass of the corresponding solution in a fixed AdS 
background (for fixed Qm, A). This inequality is of course in accord with our intuition that gravity tends 
to reduce the mass. A similar property has been noticed for monopole solutions in a spontaneously broken 
gauge theory (without cosmological constant) p9| . 

Beside these fundamental gravitating solutions, EYM theory possesses also excited solutions not pre- 
sented in a fixed AdS background. 

In this case, the metric functions of this EYM solution are considerably smaller at the origin, and 
the gauge field functions have their peaks and nodes shifted inwards, as compared to the first branch 
solutions. The energy density of the matter fields has higher peaks, which are shifted inwards, compared 
to the fundamental gravitating configurations. Otherwise, many properties of the branch axially symmetric 
solutions are similar to those of their asymptotically flat counterparts. 

To see the change of the functions for a second branch solution, we exhibit in Figs. 9 a three-dimensional 
plot for a configuration with winding number n = 2, node number k = 1, magnetic charge Qm — 1-18 and 
total mass M = 1.498. Fig. 9h presents the energy density of the matter fields e, showing a pronounced 
peak along the p-axis and decreasing monotonically along the z-axis (here x = r/(l + r). Equal density 
contours presented here reveal a torus-like shape of the solutions. 

The same general behavior is obtained for two-node solutions. 

We don't address in this paper the problem of limiting solutions, which is still unclear even in the 



spherically symmetric case. Using the metric form (52) Bjoraker and Hosotani have observed that, as the 



parameter b (in ([31])) is increased, the function H(f) hits zero from above for some values of f. Also, 
when b = b c , k has a finite value, uj(rh) = p{fh) = H'(fh) = and the space ends at r = fh- There is a 
universality in the behavior of the critical solutions |7| . The meaning of the critical spacetime is yet to bc 
clarified. 

This behavior strongly contrast with the asymptotically flat case. There are no restrictions on the 
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node number k and, as k — » oo, the BK solutions tend to a configuration that is the union of two parts. 
A non-trivial part for f < 1 represents an oscillating solution, and a simple part for r > 1, represents the 
exterior of an extremal Reissner-Nordstr0m (RN) solution with mass M = 1 and charge Qm = 1 |3C|| . The 
limiting axially symmetric configuration represents the exterior of an extremal RN solution with mass n 
and charge Qm = n fUsf . 

We have found difficult it to obtain axially symmetric generalizations of the spherically symmetric 
solutions near the critical spacetime, with large errors for the functions. A different metric parametrization 
appears to be necessary. 

5 CONCLUDING REMARKS 

In this paper we have presented numerical arguments that EYM theory with a negative cosmological 
constant possesses regular static axially symmetric solutions. They generalize to higher winding number 
the known spherically symmetric solutions. 

We started by presenting arguments that SU(2)-YM theory possess solutions with nonvanishing mag- 
netic and electric charges and arbitrary winding number when AdS spacetime replaces the Minkowski 
space as the ground state of the theory. The spherically symmetric solutions we found have properties 
similar to the lower branch of their known gravitating counterparts. 

When including gravity, we have presented results suggesting the existence of axially symmetric solu- 
tions. These configurations have continuous values of mass and non-Abelian magnetic charges and present 
a branch structure. As seen from the figures, the distributions of the mass-energy density — Tj? can be 
different from those of spherical configurations (i.e. almost toroidal distributions). As a result of these 
distributions of mass-energy density, the spacetime structure of our solutions can be considerably nonspher- 
ical, strongly axisymmetric. We have noticed a somewhat different behavior of the k = 1 fundamental 
gravitating solutions as compared to higher nodes excitations. 

We have not considered the question of stability for higher winding number solutions. However, we 
expect that the nodeless, lower branch of axially symmetric solutions is stable. A rigorous proof is however 
desirable, analogous to the proof given for the spherically symmetric case. 

In Ref. H a scaling law is derived for the mass spectrum of the spherically symmetric solutions with 
respect to their nonabelian charges Qe and Qm, A and the parameter v = AwG/e 2 . The mass of monopoles 
and dyons is expressed in terms of a universal function /(QmiQe) independent on v, A and also on k. 
The monopole and dyon solutions in the lowest branch (k = 0) are essentially the solutions in the fixed 
AdS background metric and will be stable. The solutions in the higher branches (k > 0) are obtained by 
dressing monopole and dyon solutions in the fixed AdS background metric around the BK solutions in the 
asymptotically flat space. As all BK solutions are unstable, the monopole and dyon solutions in the higher 
branches are also unstable. 

We suspect the existence of a similar behavior for the axially symmetric monopole solutions discussed 
in this paper. We have found already that for n > 1 the lowest branch monopole solutions are essentially 
solutions in a pure YM theory. Here the BK solutions are replaced with the KK solutions. The axial 
symmetry will introduce a new parameter to the universal scaling function, the winding number n. 

As discussed in H, [|, the soliton solutions depend nontrivially on the value of A. A fractal structure in 
the moduli space of the solutions has been observed |?], |3l]]. New branches emerge as |A| becomes smaller 
and the shape of branches has approximate self-similarity. As A — > solutions on a branch collapse to 
a point - the BK solution, and the nodeless solutions disappear as their ADM mass vanishes. We didn't 
investigate this point in the present work, restricting ourselves to A = —0.01. However, it is very probable 
that this general picture remains valid for the axially symmetric configurations and the continuum of 
solutions becomes a discrete set as A — > 0. Here also, in this limit, the BK solutions are replaced with the 
KK generalizations. 

Actually, due to conformal invariance of the YM equations and the fact that the AdS metric is reduced 
to the flat metric by a conformal transformation, any solution on AdS background corresponds to some 
solution in the flat space. In Minkowski coordinates these are nonstationary solutions. For example, the 



14 



monopole solution fl33|) corresponds to the well-known flat-space meron solution |32j • 

Analogous solutions with higher winding number should also exist when including a Higgs or a dilaton 
field in a theory with A < 0. We conjecture the existence of axially symmetric gravitating YM black hole 
solutions with a nonvanishing cosmological constant. These would be the AdS spacetime generalizations 
of the asymptotically flat solutions discussed in |33|. The axially symmetric nodeless solutions are of 
particular interest, because they are likely to be stable against linear perturbations. For A > 0, axially 
symmetric EYM solutions generalising the spherically symmetric configurations found in Q should exist 
also. 

Finally, this is not the complete story: the investigations can be extended to the gravitating axially 
symmetric dyon solutions. This can be an important issue, since in the asymptotically flat case this 
problem has not been solved yet. For A > 0, there are no-go theorems forbidding the spherically symmetric 
dyon regular solutions (27| [[J. Also, the authors of |^| conjectured the absence of charged regular EYM 
solutions. However, as shown in [Q the BK solutions admit slowly rotating charged generalizations. The 
total angular momentum of these solutions is proportional to the non-Abelian electric charge. Therefore it 
is natural to look for (gravitating-) axially symmetric dyon solutions, which however have not been found 
so far within a nonperturbative approach. 

In a theory with a negative cosmological constant, the boundary conditions at infinity allow for spheri- 
cally symmetric dyon solutions. Starting from the spherically symmetric n = 1 solutions, we have obtained 
nongravitating axially symmetric dyon YM configurations. We suppose that these configurations will sur- 
vive when coupling to gravity. 
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Figure Captions 

Figure 1: The gauge functions Hi and the energy density e (in units 47r/e 2 ) are shown as a function of the 
radial coordinate r for the angles 9 = 0, 7r/4 and tt/2 for three nongravitating solutions with (Qm — 2.391, 
M = 0.576), (Qm = 3, M = 0.332) and (Q M = -13.284, M = 1.118). Here the winding number is n = 3 
and A = -0.01. 

Figure 2: The gauge functions Hi and the energy density e (in units 47r/e 2 ) for a monopole solution 
with n = 3, M = 5.481, Qm = —37.185 are shown as a function of the compactified coordinates z, p. 

Figure 3: Typical nongravitating spherically symmetric solutions for A = —0.01, a fixed value of the 
parameter b = 0.001 and a = 0., 0.01, 0.02. The figure for a = corresponds to a magnetic monopole. 
The energy density e(r) is given in units 47r/e 2 . 

Figure 4: The gauge functions Hi and the energy density e (in units 47r/e 2 ) for a nongravitating dyon 
solution with n = 2, k = 1, E = 1.046, Qm — —4.124, Qe = —4.36, are shown as a function of the 
compactified coordinates z, p. 

Figure 5: The energy density e (in units 47r/e 2 ) for a dyon solution with E = 1.116, Qm = 
—2.272, Qe = 6.43 and for a monopole solution with E = 0.607 and the same value of Qm is shown 
for several values of the angle 9 (n = 2). 

Figure 6(a): The coordinate transformation between the isotropic coordinate r and the Schwarzschild- 
like coordinate f is shown for spherically symmetric solutions with k = 1 — 3. 

Figure 6(b): The metric functions /, I and the gauge field function w are shown as a function of the 
the isotropic coordinate r for the solutions presented in Fig. 6a. 

Figure 7: The metric functions /, I, m, the gauge functions Hi and the mass density are shown as 
a function of the radial coordinate r for the angles 9 = 0, tt/4 and 7r/2. Here n = 1, 2, and 3, k = 1, 
Q M /n = -8, M(n = 1) = 0.955, M(n = 2) = 2.351 and M(n = 3) = 4.179. 

Figure 8: Mass M is plotted as a function of magnetic charge Qm for first branch gravitating monopole 
solutions at A = —0.01. The winding number n is also marked. 

Figure 9: The metric functions /, I, m, the gauge functions Hi and the mass density e for a monopole 
solution with n = 2, M = 1.498, Qm = 1-18, are shown as a function of the compactified coordinates z, p. 



17 





\Q..= 2 391 


i 




/ \ 
/ x 

/ \ 
/ \ 


9= 7C/2 




/ \ 
/ \ 
i , 
/ x 
/ \ 


e= 7c/4 




/ ~ \ 


0=0 




/ Q M =3 \\ 

' t \ \ 
/ ' ^\ \ 






Q M ="1 3.284 











12 3 4 

iog 10 (r) 

Figure la. 



18 




19 





1 


1 1 

\ Q M = -1 3.284 




/ 

/ 

/ 

/ 

/ 

/ 


e= 7c/2 




/ 

/ 

/ 

/ 

/ 

/ 


e= 7c/4 




/ 
/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


e= o 




/ 

/ 

/ / 


Q M = 2.391 \ 


i i i 



12 3 4 

iog 10 (r) 

Figure lc 



20 




12 3 4 

iog 10 (r) 

Figure Id 



21 




22 







Figure 2a. 



23 




Figure 2b. 



24 




Figure 2c 



25 




Figure 2d 



2G 



£ 

OjOODG r 



OjOOOi - 




Figure 2e 




27 



0.5 







-0.5 



1 — ■ — 


i 

■ v 




i i 










e(r)x1000 




\ \ 




E=0.035, Q M =0.9, Q E =0 








E=0.43, Q M =0.58, Q E =-2.26 
E=1.11, Q M =0.69, Q E =-3.31 

- • — "' 




\ \ ' 




i 




u(r)-1 



-10 12 3 4 

iog 10 (r) 

Figure 3 



28 




29 




Figure 4b 



30 




Figure 4c 



31 




Figure 4d 



32 



H 5 




Figure 4e 



33 




Figure 4f 



34 



£ 

O.OOl r 



00005 - 




Figure 4g. 



35 



1 

log-io(r) 



Figure 5 



36 




-2-10 1 2 

log 10 (r) 

Figure 6a 



37 




-2-1012 

log 10 (r) 

Figure 6b 



38 



1 


1 


i i 

e= 7c/2 




n=3 \ 


e=7u/4 




i \ 
j \ 


e=o 




i i \ \ 

n= 2 

i 





iog 10 (r) 



Figure 7a 



39 




40 




Iog 10 (r) 



Figure 7c 



41 




42 



0.9 



0.8 



1 


/ / . ' /■' i 
/ / ' / '' ' 


s= s T-"-~ 1 


n= 1 




e=7i/2 




/ ' 1 1 ' 
/ ' 1 1 ' i 
/ • . i 
/ ' f 1 1 i 

/ / i ' i 1 

/ / 1: i 
/ e ' 


e=7i/4 


n=2 


1 ' 
1 1 1 


e=o 


n=3 y 


/ / ' 

/ / i 
/ i 




i 


i 





1 2 

'og-io(r) 

Figure 7e 



43 




44 



0.95 



E 

0.9 



0.85 





1 

n= 1 


/ / f! ' ' 

/ / il i ■ ■' 
/ / / 1 ' 1 

// '' 1 


^^^-1 \ 

e= 7c/2 






/ 'I 

/ // I " ' 


e=7c/4 




_n= 2 


J ' ' 


e=o 




n= 3 


/ / 1 






i 


1 





1 2 

log 10 (r) 

Figure 7g 



45 




46 




47 



hi 




Figure 9a 



48 




Figure 9b 



49 




Figure 9c 



50 




Figure 9d 



51 




Figure 9e 



•52 




Figure 9f 



53 




Figure 9g 



54 




55 



